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Recent results in light hadron spectroscopy are reviewed. Attention is given to the requirements of precision 
determinations in lattice gauge theory. Different methods for extracting the running coupling constant as are 
compared. Some recent results in determining the nature of the interaction between static quarks are also 
discussed. 



1. INTRODUCTION 

In principle lattice Monte Carlo methods are 
capable of precision evaluations of the QCD spec- 
trum. Furthermore one can study the effect of 
varying the fundamental parameters, namely the 
quark masses, from their experimental values. 

In order to organise the discussion, I present 
the current state of QCD spectrum calculations 
by comparing them to the ideal target of perfect 
accuracy. The topics to be covered are the control 
of the discretisation errors, finite volume errors, 
quark mass errors and ground state extraction 
errors. 

In a later section I discuss the link to pertur- 
bation theory. This is relevant for an accurate 
extraction of the running coupling of QCD. It is 
also relevant for other perturbative contributions, 
such as the matching coefficients (Z) needed to re- 
late lattice matrix elements to continuum matrix 
elements. 

The interaction between static colour sources 
is discussed since recent results shed light on the 
confinement mechanism. 

2. PRECISION QCD ON A LATTICE 

One aim of lattice QCD is to be able to repro- 
duce accurately the hadron spectrum. This will 
then give confidence in the method and will allow 
evaluations of new quantities of interest such as 
matrix elements. A study of the effect of varying 



*To be published in the Proceedings of Lattice 1994, NucL 
Phys B (Proc. Suppl.), Liverpool Preprint LTH 340, hep- 
lat/9412032, December 5, 1994 



the quark masses is also of great interest. Because 
of the extreme computational resource needed for 
full QCD simulation, progress can be made by 
studying quenched and partially quenched sys- 
tems as well. 

There have been new dynamical fermion calcu- 
lations reported this year. A growing trend 
is to compromise between quenched and dynami- 
cal calculations by using partial quenching. This 
implies using valence quarks with different prop- 
erties (mass, or discretisation) from the dynam- 
ical (sea) quarks present in the vacuum. The 
quenched approximation itself is still a rich sub- 
ject of study and much remains to be done. 

The ultimate aim of lattice QCD calculations is 
to extract physical results applicable to the con- 
tinuum limit. To organise the presentation of the 
current state of the subject, I consider in turn the 
various approximations inherent in the lattice ap- 
proach. 

2.1. Discretisation Errors 

The lattice spacing a is introduced to regulate 
the field theory of QCD. Then dimensionless ra- 
tios of mass values determined from lattice cal- 
culations will have errors of order Ma or Aa, in 
general, where M is a mass value and A represents 
a typical QCD momentum scale. The continuum 
limit is obtained by extrapolating such ratios to 
a — 0. Even if only states with light quarks (u,d 
and s) are studied, the masses can be large: the 
Sl~ has a mass of 1.67 GeV. So, a priori, correc- 
tions of order Ma can be large for current values 
of the lattice spacing of a^^ < 3 GeV. 
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Note that dimensionless ratios such as M/A 
will have discretisation errors of order or log(a) 
since A is defined through matching with pertur- 
bation theory. Thus one should not use such ra- 
tios to obtain a continuum limit. 

For the study of hadrons made of quarks, the 
Wilson discretisation of fermions is simple and 
flexible. It has errors of 0{a). As an illus- 
tration, consider the nucleon to p mass ratio. 
The most extensive results have been obtained 
in the quenched approximation. Then keeping 
the lattice volume constant in physical units (at 
MpL K, 9), and after extrapolating the quark 
masses to the physical light quark masses, the 
ratio has been fitted [01 as 



N 



1.28(7) 



0.27Mpa 



(1) 



as illustrated in fig This rate of a-dependence 
implies that an extrapolation in a is needed even 
with data at /3 = 6.17 on 32^ x 30 x 40 lattices 
where MpO w 0.3. The results shown in fig. 1 
from the GFll group ||^ are the most compre- 
hensive available. Nevertheless, the statistical er- 
rors are quite large so that, even in the quenched 
approximation, results for Mj^/Mp are not yet 
precise. Thus although the GFll data are con- 
sistent after extrapolation in a and L with the 
experimental ratio of 1.22, much work still needs 
to be done to get an accurate calibration of the 
quenched approximation value. A recent analysis 
by the LANL group Q using 32^ x 64 lattices at 
/? = 6.0 has also produced relevant data and this 
is shown in fig 1. Even with 80 configurations, 
they are unable to reduce their errors sufficiently 
to provide a stringent cross check. 

The situation for staggered fermion discreti- 
sation in the quenched approximation is sum- 
marised in the review by Gottlieb g. A gam 
current data are consistent with an a-dependence 
which yields a continuum limit close to the exper- 
imental ratio of 1.22. 

Given the substantial 0{a) corrections ob- 
served with Wilson fermions, much effort has 
gone into formulations that improve the discreti- 
sation errors. The SW-Clover fermionic action 
has 0{asa) errors with a correct treatment of ob- 
servables Simulations using this action |^J^ 
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Figure 1. The nucleon to p mass ratio from 
quenched results extrapolated to the chiral 
limit. 



at P — 6.0 and 6.2 are not yet of sufficiently high 
statistics to enable Mjsi /Mp to be determined ac- 
curately. In fact I return to discuss the current 
state of the UKQCD nucleon data when I discuss 
ground state extraction later. 

For heavy quarks, the hadron spectrum is even 
heavier and 0{Ma) corrections are even more im- 
portant. For c quarks, the SW-clover formula- 
tion is promising and has been successfully used. 
For b quarks, it is feasible to use a formalism 
which treats the heavy quark explicitly: non- 
relativistically (NRQCD), as static quarks, or by 
a discretisation of heavy quark effective theory. 
Such heavy quark calculations are covered in talks 
by Sommer [||, Sloan Q etc. 

Another approach to improving the discretisa- 
tion errors, is the mean-field improved formal- 
ism fll|| . This provides an impressive description 
of many features. What is needed is a method to 
estimate reliably the residual error. 

The goal of removing completely the discreti- 
sation errors is achieved by the perfect action. 
For the gluonic sector of the theory, the standard 
Wilson plaquette discretisation already has errors 
only 0{a^). Current approximations to the per- 
fect action are reported by Hasenfratz [n2|. As 
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well as a perfect action, perfect observables will 
be needed to remove lattice artifacts such as the 
lack of rotation invariance in the static potential 
for example. 

2.2. Finite Volume Errors 

Corrections in a spatial volume can arise 
because of exchange of the lightest hadron (the tt 
meson of mass m^) around the boundary. Finite 
size corrections can also occur if the radius of a 
hadronic state is not small compared to L. To 
remove these corrections, one needs 



TO. 



rL » 1, and rh/L « 1 



(2) 



The finite value of L also makes momentum dis- 
crete and for some studies it is appropriate to re- 
quire the minimum momentum qmin ~ 271 / L « 

TO^. 

There can be surprises: in pure gauge stud- 
ies the finite size effects on the glueball spectrum 
are significant out to mcL « 8 where ttiq is the 
lightest state (the 0++ gluebaU). So ">>" really 
means ">>" not just ">" in this case. This is 
understood and arises since the dominant finite 
size effect comes from flux loops encircling the 
spatial boundary (called spatial Polyakov loops 
or torelons). 

' ' " to 



Moreover, these effects were predicted |13 



be approximately twice as strong in full dynam- 
ical quark simulations since one torelon (rather 
than two) can then mix with the glueballs. 

At the Lattice 1993 conference, results were 
presented showing that finite size effects on 
the hadron spectrum indeed appear worse for full 
QCD compared to the quenched approximation. 
The origin of the effect can be seen from the ex- 
ample of a meson correlation illustrated in fig |^ 
where one quark propagator encircles the spatial 
boundary. So effectively, the quark feels the (non- 
abelian) phase of the spatial Polyakov loop. Now 
in the quenched approximation, the symme- 
try ensures that such loops have zero expectation 
value. In contrast, the quark mass term in the 
dynamical quark case breaks this symmetry. 
This then allows the Polyakov loops to become 
non-zero on average and so enhances their contri- 
bution to finite size effects. 




Figure 2. A finite size effect 



A new study of this topic was reported by 
the APE collaboration ||]. They use anti- 
periodic spatial boundary conditions on the Wil- 
son dynamical (sea) quarks. This 'polarises' 
the Polyakov loops so reducing statistical errors. 
They then explore the magnitude of the resid- 
ual finite size effects by measuring tt and p meson 
masses using both periodic and anti-periodic spa- 
tial boundary conditions for the valence quarks. 
With L « 1.5 fm, they find finite size effects on 
the masses of only a few per cent. 

2.3. Quark mass extrapolation 

The requirement is to incorporate fully quarks 
of physical mass in the QCD vacuum. Because 
of the large computational resource needed, typ- 
ical studies are limited to spatial sizes of 16^ and 
quark masses corresponding to mj^/mp w 0.3. 
For technical reasons, one of the preferred meth- 
ods is to include two flavours of degenerate light 
quarks using the staggered fermion discretisation. 

In most simulations, the main effect of in- 
cluding dynamical quarks is a shift of the /? 
value. Thus, instead of comparing at the same 
bare (3 value, one should compare at the same 
lattice spacing (ie same ma value). Then ra- 
tios of masses appear very little changed com- 
pared to the quenched case. Since the present 
quenched spectrum is consistent with the exper- 
imental data, it is to be expected that adding 
dynamical quarks will not change this state of af- 
fairs. One expected effect of dynamical quarks is 
to modify the small distance forces. This region 
is amenable to perturbative evaluation and, be- 
cause of the factor 33 — 2Nf in the denominator 
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of the running coupling, a stronger interaction is 
expected at short distance. This can be tested 
by looking at the force between static quarks at 
short distance |15|. 

An extrapolation of the light quark masses to 
the physical value is needed. For full QCD, data 
are not sufficiently precise to explore this extrap- 
olation accurately. Further work is needed to 
check whether the extrapolation is smooth: some 
"threshold effect" might occur as the quark mass 
is reduced (the opening up of the decay p ^ tt + tt 
has been suggested). 

For studies of the quenched (and partially- 
quenched where the valence quarks are different 
in mass from the dynamical or sea quarks) ap- 
proximation, the limit of small quark masses is 
potentially tricky. This has been studied in gen- 
eralised chiral effective theories and is reported 
on by Gupta |16|. To circumvent these prob- 



lems associated with extrapolation to the chiral 
limit, one way is to present lattice results at the 
quark mass used. This is appropriate for com- 
paring between lattice calculations (eg. the Ed- 
inburgh plot) but may also be used to compare 
with experiment in principle. From lattice results 
for hadron masses for a range of quark masses, it 
should be possible to develop an accurate semi- 
empirical mass formula relating the hadron mass 
to the masses of its quark ingredients. Then ex- 
perimental data for the combinations of quark 
masses (u, d and s) accessible can be used to fix 
the coefficients in such a model. Thus lattice data 
would be used to validate a quark-model mass 
formula which could then be used for direct com- 
parison with lattice results. 

While computational resources limit full QCD 
work, there is still a lot of physics that can be 
studied in the quenched approximation. The 
quenched approximation has spontaneous chiral 
symmetry breaking and gives a surprisingly good 
overall description of the QCD spectrum. It is 
even possible to estimate matrix elements within 
the quenched approximation which enable one to 
study processes beyond that approximation. Ex- 
amples are the decay p — > tt + tt; glueball decay 
and the 77' mass. 

Here I discuss some of these analyses to illus- 
trate the possibilities. 




Figures. The glueball tttt correlation 



2.3.1. Glueball decay 

The glueball spectrum has been studied in the 
quenched approximation p^ , p^ and the lightest 
glueball is found to be the — 0+"'" state at 
around 1600 MeV with the next heaviest state of 
2++. 

Within the quenched approximation, the light- 
est glueballs cannot decay. In full QCD, they 
can decay to mesons: for instance a pair of pseu- 
doscalar mesons. There will also be mixing be- 
tween glueballs and mesons with quark content. 

There have been 20 years of unsuccessful exper- 
imental attempts to identify a glueball. To guide 
this search, it is important to estimate the mag- 
nitude of the matrix elements between glueballs 
and light quark mesons. Then, should the ma- 
trix elements be small, the assignment of glueball 
status could be made with confidence. Candi- 
dates for the 0+"'' glueball are the /o(1590) meson 
with width 180±17 MeV and the /j(1710) meson 
(was called 0(1690)) which has a decay width of 
140 ± 12 MeV. 

From the point of view of lattice evaluation, 
such matrix elements are hard to measure because 
of a mis-match between the techniques used to 
extract glueballs and those used for light quark 
mesons. A clear glueball signal can only be ob- 
tained by analysing many configurations (typ- 
ically thousands of independent configurations 
with all time-slices used as a source). Because 
of the computational overhead of inverting the 
fermionic matrix, quark propagators are usually 
only evaluated for hundreds of configurations and 
from only one source per configuration. Further- 
more the glueball signal only extends for a few 
time steps whereas a pseudoscalar meson can be 
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Figure 4. The glueball decay coupling. 



followed out to very large time. 

A pioneering attempt to tackle this problem 
has been made by the GFll group They 
use lattices of spatial size 16^ at /3 = 5.7. Us- 
ing Wilson quarks with K = 0.165, they measure 
3800 configurations. They are able to isolate the 
O"*""*" glueball, pion and two pion states. The two 
pion state is calibrated by exploring expectation 
values of 4 pion operators. Then the matrix el- 
ement between a glueball and such a two pion 
state is studied - see fig |[ In principle the inter- 
mediate state has many different contributions: 
glueball, excited glueballs, 7r(A:) + 7r(— /c) with rel- 
ative momentum k = 2TTn/L and 0^"*" mesons. 
To separate these contributions by studying the 
t-dependence of the observable is difficult in prin- 
ciple and hopeless in practice. Since a reasonable 
signal only survives out to t = I with their cur- 
rent statistics, the GFll group can make, at best, 
a very rough estimate of the size of the matrix ele- 
ment. Even such a rough estimate can be useful: 
for example by indicating that a large glueball 
width is not allowed. In the quenched approxima- 
tion the glueball does not decay and, moreover, 
the analysis is performed with a choice of quark 
mass such that m{0~^'^) « 2mT^. Thus the on-shell 
glueball decay matrix element can be estimated 
from the lattice measurements. 

The result shown in fig is consistent, within 
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Figure 5. The ratio of disconnected to connected 
diagrams. 



large errors, with that from the observed decay 
of the /j(1710) meson. Some evidence is also 
obtained that the decay matrix element increases 
with increasing pseudoscalar mass - as observed 
experimentally. See ref jl^ for more details. 

2.3.2. The r/' mass 

In the quenched approximation, there are no 
quark loops in the vacuum. Then, since the vac- 
uum has isospin zero, the flavour octet and singlet 
pseudoscalar mesons (rys and 771) will be degener- 
ate if all quark masses are the same. Taking the 
strange quark mass into account from the tt and 
K leads to an rjs at 567 MeV and an r/i at 412 
MeV. Experimentally the mixing is such that ap- 
proximately the ?7g is the 77(547), while the rji is 
the heavier 77' (958) meson. Thus in the quenched 
approximation, the contribution from the rj' will 
have a mass which is much too low. This, in- 
cidentally, is at the root of the modified chiral 
behaviour in the quenched approximation. 

Within the quenched approximation one may 
study the mixing between quark loops and the rj 
mesons explicitly. Let Pg — {p^ + '^tg)"^ be the 
quenched 77 propagator (which is the same as the 
77' propagator in the quenched approximation), 
Then consider the effect of a quark loop Q on the 
77' propagator Pi = (p^ + ttij)^^: 

Pi = Pg - PgQPg + PgQPgQPg + . . . 



Ps^ +Q + ml + ml ^ ^ 

Thus — rug + m'^ where uiq is the matrix ele- 
ment of a quark loop Q in the quenched vacuum. 
This mixing works to increase the mass of the 
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Figure 6. The nucleon effective mass versus t 
from UKQCD LL data witfi one and two expo- 
nential fits. 



Figure 7. The nucleon effective mass versus t 
from UKQCD preliminary data for LL (squares), 
SL(crosses) and SS(diamonds) with a two expo- 
nential fit. 



singlet - so giving a chance to obtain the experi- 
mental rj' — r] mass splitting. Lattice methods to 
explore this have been proposed The basic 
idea is to compare the disconnected r]'{0)r]'{t) cor- 
relator with the connected correlator - see fig |^. 
The ratio is proportional to the required mixing 
matrix element mp. The technical difficulty is 
to measure the disconnected correlator: an in- 
genious method was used to cut down the num- 
ber of propagator inversions by putting sources 
at each spatial site on a timeslice and allowing 
gauge invariance to select the required contribu- 
tion. An analysis using Wilson valence quarks at 
/3 = 5.7 on 12^20 and 16^20 lattices yields mix- 
ing Too = 751 ± 30 MeV to be compared to the 
experimental value 

Too = ml, -ml^ (852 MeV)^ (4) 

So indeed quenched QCD does yield a reasonable 
understanding of the 77' mass. 

The Witten-Veneziano formula relates the ij 
mass mixing to the topological susceptibility x 

by 



Measurements of the topological susceptibility in 
lattice gauge theory are plagued with uncertain- 
ties from lattice artifacts. Using the cooling 
method on the same lattices as above, a result of 
ml = (1146± 67MeV)2 was obtained [|o). This 
is qualitatively similar to the value obtained by 
studying the mixing directly and is further con- 
firmation that the underlying physics of the rj' is 
understood. 

2.4. Ground state extraction 

The euclidean time formulation implies that 
eigenstates of the transfer matrix contribute to 
hadronic correlators <H'^{0)H^{t)> as 

Cat{t) = <H-{0)HHt)> 

hlh\{e-"'^^ + e-™i(-^-*)) + . . . (6) 

where hi is the amplitude to produce eigenstate 
i from operator H'^. The expression is writ- 
ten assuming periodic boundary conditions in 
time for bosonic eigenstates. The data is usually 
presented by computing an effective mass from 
"loff (i) = - \og{C{t)/C(t - 1). The ground state 
contribution dominates moff(t) at large t since 
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mo < rrii for all i 7^ 0. In practice, extracting 
the ground state mass mo from data with statis- 
tical errors which increase with t is subtle. I wish 
to emphasize that a good procedure is: 

(i) a multi-exponential fit to the widest accept- 
able t-range with 

(ii) several hadronic operators to stabilise 
the fit. 

One way to understand this guide is that the 
ground state is only determined accurately when 
an estimate of the first excited state is available. 
This is necessary since the energy difference con- 
trols the rate of approach of Cabit) to the ex- 
pression given by the ground state component 
alone. However, fitting 2 (or more) exponentials 
to just one function C[t) is not very stable: better 
is to have several such functions (provided that 
they do indeed have different relative amounts of 
ground state and excited state). 

As an example of this discussion, I use some 
preliminary UKQCD data for the nucleon mass 
determination. With just one quantity measured 
(the local source - local sink: LL), a plateau in 
the effective mass was found as shown in fig ^ 
This shows a typical problem: the effective mass 
decreases out to large f-values. I have called this 
the "droop". It is difficult to be sure whether 
such a behaviour is statistically significant or not. 
For instance, a 2-exponential fit to the same data 
gives a significantly higher value of mg. A reso- 
lution of this dichotomy comes from adding more 
measurements: SL and SS where S refers to a 
smeared source or sink. The combined fit to the 
3 types of data is shown in fig |^ and strongly 
favours the 2-exponential fit with the higher mg 
value. 

2.4.1. Correlations among data 

It is widely known that, particularly for 
hadronic measurements, the data Cabit) are 
strongly correlated statistically. For example 
Cab{t) and Cab{t + 1) are often 99% correlated. 
By this I mean that each fluctuates among lattice 
samples within its error, but that when the value 
at t is above average then so will be the value 
at i -I- 1 . This means that each t value is very far 
from statistically independent. The standard way 
to accommodate this is to use 'correlated x^'- 



The drawback of this method is that it can be- 
come unstable if the number of lattice samples N 
is insufficient when there are D types of data be- 
ing fitted (ie D is number of t- values times num- 
ber of types of measurement). The expected 
per degree of freedom from such a correlated 
fit is 

^^l + ^+0(A^-) (7) 

Moreover, ii D > N, then becomes infinite 
since the correlation matrix among the data is 
non-invertible. 

This is a problem in practice: if 15 i-values are 
to be fitted with 3 types of observable then D = 
45 while N is typically between 50 and 200. The 
recommended way forward is to parametrise the 
correlations among the data with fewer parame- 
ters than the general case (which has D{D—l)/2). 
This should ensure a stable evaluation of the in- 
verse of the correlation matrix and, in particular, 
avoid the problem of spurious small eigenvalues 
of the sample correlation matrix. Various sugges- 
tions have been made to accomplish this: an SVD 
inverse 3 and 5-diagonal inverses aver- 
agingof small eigenvalues of the correlation ma- 
trix Ig^, etc. The methods used will need to take 
into account the nature of the correlations among 
the data in each application and then model them 
in a stable way p3| . A program incorporating a 
broad selection of these methods is available from 
the author. 

The alternative to such modelling of correla- 
tions is to have enough independent data samples 
N. Data sets containing of the order of thousands 
of configurations analysed are going to be needed 
for precision hadron spectroscopy on a lattice. 

2.4.2. Efficient hadronic operators 

The strategy for fitting the data to extract mass 
values involves using several hadronic operators. 
For example local and smeared operators were 
used in the previous illustration. The availabil- 
ity of several operators allows a variationally im- 
proved operator to be constructed for the ground 
state. It is, never the less, necessary to select such 
an operator basis wisely. 

It is important to explore this choice of oper- 
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Figure 8. The Bethe Saltpeter wave function for 
the p and p' from ref 



ators to find efficient ones. The definition of ef- 
ficient in this context is that the contribution of 
excited states is minimised relative to the ground 
state contribution. Non-local operators modelled 
on heavy quark bound states are promising can- 
didates: they would be the optimum choice in the 
limit of heavy quarks and so may still be appro- 
priate for lighter quarks. The construction used 
is to join the quark and anti-quark in a meson by 
a straight colour flux path of length R. Following 
experience of potential studies (24j , this flux path 
is made from smeared links. The resulting overlap 
amplitude < 0|g(0)f7(0 ^ R)TMq{R)\M > is of- 
ten referred to as the Bethe-Saltpeter amplitude 
of meson M (where Tm is Ji for the p, etc). This 
has been studied before |^,^ to learn about the 
spatial distribution of hadrons. The new obser- 
vation is that the same quantity for the first 
excited state (eg. p') has a node. Thus choosing 
the separation R at this node will give an effi- 
cient operator to make the ground state with no 
contamination from that first excited state. As 
illustrated in fig ^, this node is at i? « 3GeV~^. 
Results for baryons are also qualitatively simi- 
lar [§1. 



The world average of the QCD running cou- 
pling compiled by Webber is 
7 0-117(5). 

This average includes input from lattice analyses. 
Indeed the lattice determinations have the small- 
est errors quoted. Since the precise value of as 
is important for phenomenology and for attempts 
to look for evidence of models beyond the stan- 
dard model, it is important for the lattice com- 
munity to understand fully the source of errors 
in the lattice determinations. This is highlighted 
by the change from the value 0.108(6) quoted 
at Lattice 1993 to a value of 0.115(2) quoted | 30|] 
this year. 

There are two requirements to measure as on 
a lattice: 

(i) an accurate match to perturbation theory 

(ii) an accurate lattice energy scale set by some 



experimental input (eg. 
onia, etc) 



V{R), IP-IS heavy 



3.1. Perturbation theory on a lattice 

First I discuss the requirement that perturba- 
tion theory is accurately reproduced on the lattice 
so that as can be defined. Consider first pertur- 
bation theory in terms of the bare lattice coupling 

= 6/(3 with a = g'^/An. The plaquette action 
Sa has been calculated for iV^ = to order 

for the Wilson gauge action. For SU(3) colour, 



an 



--log(5o) 
a(l + 3.37q; -( 



17.69^2 + ...) 



(8) 



Now at /3 = 6.5, a — 0.07346 and the last ex- 
pression evaluates to 92% of the first expression 
obtained from the measured lattice plaquette ac- 
tion. The a^ term contributes 7%. So we do not 
have a good convergence of perturbation theory. 
Increasing P substantially to decrease a is not 
feasible with present computing techniques. 

The way forward has been known for a long 
time. A physical definition of a should yield a per- 
turbative series with coefficients which are much 
smaller than in the case of the bare coupling. 
Consider several such physical definitions: 

(i) av from V{q) or Sn [0. 
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(ii) aqq from the force dV/dR |3^ . 

(iii) asF from the Schrodinger functional jss). 
A variant of asp with twisted boundary condi- 
tions has recently been used jsij . 

Consider the first case, where the Fourier trans- 
form of the static potential is used as a definition: 
?>q^V{q) = —IQiravisi)- Then the logarithm of 
the plaquette action can be re-expressed in terms 
of this definition: 

an = av{l - (1.19 + Q.Q7Nf)av + 0{al)) (9) 

where ay has scale q = 3.41/a Indeed the 

coefficient of is now 0(1). The expectation is 
that all coefficients in such an expansion are 0(1). 
For our present purposes it is the coefficient of the 
ay term above that is assumed to be 0(1). 

This coefficient can be obtained by evaluating 
the perturbative connection between the bare and 
y— scheme coupling to order since the link be- 
tween Su and bare coupling is known (though 
only for Nf = 0). Indeed, since the calculation 
of the connection between and the bare cou- 
pling is under way |Q , it will soon be possible to 
determine a-^ directly from Su to relative accu- 
racy of O(a^). Meanwhile, the connection of the 
MS and V schemes is only known to lower order. 
For example for SU(3) of colour, 

a,^ = ay(l-(0.822-0.0887V/)av-|-O(a^))(10) 

When comparing a in different schemes, the 
scales in the two schemes can in principle be cho- 
sen differently. Possible choices are to take them 
equal, choose them so that the O(a^) term is ex- 
actly zero, or some compromise. From a compar- 
ison of the Schrodinger Functional and Twisted 
Polyakov methods in SU(2), there is a definite 
preference for the prescription of choosing 
the relative scale (as the ratio of A parameters) 
so that the the 0{a^) term is exactly zero. As a 
further example of this ambiguity, eq |l^ can be 
expressed Q as 

"ms(9) = av{e^'^q){l + 0.637ay + 0(a^)) (11) 

The difference between these expressions is O(a^) 
and so a preference can only be made when the 
calculation is completed. 

Until the calculations are completed, some 
cross checks are worth making. One approach 



is to use Monte Carlo evaluation on lattices with a 
wide range of (i to estimate the size of the ay cor- 
rection in eq. (^ . The method is to compare the 
extraction of ay from 1x1, 1x2, 1x3 and 2x2 
Wilson loops. Then, relative to the (unknown) 
coefficient of the 1x1 extraction, the other cases 
have a^ coefficients of 0.36, 0.53 and 0.97 respec- 
tively. This is indeed confirmation that coefh- 
cients 0(1) are plausible. 

Another cross check is to use the recent pertur- 
bative evaluation [3^ ] of as f in terms of the bare 
coupling to obtain for SU(2) of colour: 

an = asF(l - 0.737asF + 0.039a|^ + . . .) (12) 

with scale a — q^^ = L. This is the first rela- 
tionship to be obtained to this order between two 
physically defined lattice couplings. The coeffi- 
cients are indeed less than one, although with the 
caveat that SU(2) colour tends to have smaller 
coefficients than SU(3). 

Besides evaluating ay from Su, it can be 
directly obtained from the static quark poten- 
tial IHJ. Present results are compatible between 
the two methods within the inherent errors of 
each. 

Many phenomenological models have been pro- 
posed to obtain an improved definition of a in 
lattice studies ||40|-|4^] by modifying the bare cou- 
pling. These suggestions were constructed with 
the intention of rendering the perturbative coeffi- 
cients of 0(1). In this respect they are similar to 
the av scheme outlined above but, because they 
are empirical in origin, they are less appropriate 
to a precision analysis of a. 

Since Su is the shortest distance quantity de- 
fined on a lattice, it corresponds to the highest 
momentum. It is also easy to measure accurately. 
Thus it is a very popular starting point for ex- 
tracting a. Indeed, a-^ can be obtained from 
the measured Su with error 0{a^). In order to 
check that perturbation theory is accurate, the 
coefficient of this a^ term should be evaluated 
(this is in progress) and then the values of a-^ 
from different lattice spacings should be com- 
pared with the perturbative evolution through 
the beta-function. If this perturbative evolution 
were accurate above some scale q = 1/a, then 
one would have confidence in this method of ex- 
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tracting a. For the other methods of defining 
and extracting a, the same checks are needed 
and they have been thoroughly explored for the 
Schrodinger Functional scheme in the quenched 
approximation ]3^ . This check of the evolution 
is essential to control possible non-perturbative 
effects which I now discuss. 

3.2. Non-perturbative contributions 

Since the method employed is to extract a 
from lattice measurements, there may be a 
non-perturbative component present in princi- 
ple. Note that such non-perturbative com- 
ponents are also equally important in contin- 
uum/experimental determinations of 0:5. I shall 
illustrate this problem by considering the "force" 
definition of the coupling. For SU(3) colour with 
Nf = 0, this is 



dV 



1 



R 



(13) 



Now, since the force dV/dR has a non- 
perturbative contribution from the string tension 
if, one could equally well introduce a modified 
definition by subtracting it: a = aqq — 3KR^ /A. 
Since a « bo/\og{R~^), the two definitions of a 
differ by a term proportional to exp(— 60 /a) . This 
is a non-perturbative difference: the two defini- 
tions agree to all orders of perturbation theory. 
For q — R^^ >> 10 GeV, this non-perturbative 
piece is negligible (< 1%). However, one should 
keep in mind that different "physical" definitions 
of a could have different non-perturbative compo- 
nents. For example when extracting ay from Sn, 
the typical momentum scale is large (q* ~ 3.4/a) 
but it arises from a cancellation of = and 
qi = ±7r/a. So non-perturbative effects could 
arise even though q* > 10 GeV. 

3.3. Comparison of as determinations 

I now make a comparison of several methods of 
extracting a-^ from lattice measurements. To 
make the comparison as close as possible, the 
scale is set in each case from the best measured 
quantity (the static quark potential at moder- 
ate interquark separations) with the definition 



r'^dV/dR\ro = 1-65. This definition |43| corre- 
sponds to ro « 0.5 fm and to Vli'ro « 1.2. To 



have the most accurate lattice data, I use the 
quenched determinations for SU(3) colour. 

(i) Then at q = 37/ro (approx 15 GeV), the 
Schrodinger functional method obtains 

a^g= 0.1108(23)(10) 
where the first error is from the scale and the sec- 
ond is assuming the correction has coefficient 
±1. 

(ii) The method of extracting the coupling from 
the plaquette action 5*0 (via ay) yields 

ajjrs = 0.1151(5)(13) 
where data |46| at j3 = 6.5 has been used which 
corresponds to q* = 38.3/ro and this has been 
run to 37/ro using 2-loop formulae. 

(iii) The last method uses the force and, at the 
highest /3- values currently available of 6.5, 
corresponds to q < 7/ro only. However, it gives 
results for a range of q-values consistent with 
two-loop running. Assuming two-loop running to 
37/ro, then yields 

a- =0.1180(21) (13)(+g°) 
where the last error term comes from the possi- 
ble non-perturbative effect of removing the string 
tension before running the coupling. There is no 
sign in the data for any such non-perturbative 
effect but its absence can only be confirmed by 
going to smaller lattice spacing. 

These values are in reasonable agreement with 
each other. The 0{a^) terms need to be calcu- 
lated to refine the comparison further (this will 
remove the 1% error from assuming the coeffi- 
cients are within ±1). A guide to the possible 
size of the 0(0;'^) terms comes from comparing 
the different definitions of a at different scales 
so that the terms vanish - rather than at the 
same scale. Some evidence exists [ p^ , ^8|j35| that 
this prescription works well for SU(2). For the 
Schrodinger functional method, this would then 
yield 

aj^g = 0.1141(24)(15) 
The difference from the result quoted above is ap- 
proximately 2.5a^ which shows that coefficients of 
that magnitude are plausible in this case. More- 
over, this scale choice gives a value in closer agree- 
ment with the other methods. 

After this perturbative calculation, any re- 
maining difference in a values must be ascribed 
to non-perturbative effects. Since any non- 
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perturbative effects are rather different in tlie first 
method (which involves the small-volume vac- 
uum) than in the latter two determinations, this 
will be a valuable cross check. 

This comparison highlights that caution is 
needed in extracting a from lattice data. Con- 
sistency as the momentum scale is increased is 
essential. 

3.4. Setting the Scale 

The requirement is that dimensionless ratios of 
physical quantities are independent of the lattice 
spacing a so that the scale can be set unambigu- 
ously. Because of uncertainties in corrections due 
to quenching and to 0(a) effects, there is still a 
measure of uncertainty in setting the scale. 

The most accurate lattice measurement is in- 
deed To from the static potential V{K) as dis- 
cussed above. This definition, unlike that in 
terms of the string tension, is also applicable to 
full QCD. Indeed this is the method of choice for 
comparisons between different lattice evaluations. 
Of course, ro is not directly known experimen- 
tally, so other quantities are also used. 

The p mass is often chosen although the extrap- 
olation of Tfiq — !■ (or the physical light quark 
masses) introduces errors. The effect of incom- 
plete treatment of the decay p — > tt -I- tt is most 
probably quite small on the p mass. 

Another choice is /,r but the renormalisation 
factor Za then needs to be evaluated. This 
needs an accurate perturbative evaluation - cross 
checked by non-perturbative methods - see the 
review by MartineUi 

A very popular choice is the IP- IS mass split- 
ting in heavy onia. This choice is supported by 
two encouraging indications: (i) empirically, the 
mass splitting depends very little on the heavy 
quark mass and (ii) the relevant momenta (for 
instance as indicated by the wave function) are 
expected not to be so sensitive to very light 
quark effects. The drawback is that for relativis- 
tic fermions m,ca is comparable to 1, so 0{a) ef- 
fects are potentially large. This can be avoided 
by using an effective lagrangian appropriate to 
heavy quarks for the hh system: for example a 
non-relativistic treatment. This requires a care- 
ful calibration of the effective lagrangian. The 



NRQCD method, see ref [0, involves an effec- 
tive lagrangian which can be refined in principle 
by evaluating more terms. The inherent errors 
are not straightforward to estimate. 

For the quenched approximation, the (IP- 
1S)/(2S-1S) ratio for hi comes out |§ from 
NRQCD as 0.71 (similar to the static potential 
value [p4[) which differs from the experimental 
value of 0.78. This can induce some uncertainty 
in setting the scale. 

The running of the coupling reduces the depen- 
dence on the scale as the momentum is increased: 
Aa/a w 0.2 Ago/go at g = Mz- Thus a 10% er- 
ror in the scale go results in a relative error of 2% 
on a at Mz- 

3.5. as from dynamical quark simulations 

Until recently, attempts to evaluate ag from 
lattices used quenched simulation and corrected 
for the effects of quenching via a phenomenolog- 
ical potential approach |2^. This certainly cap- 
tured one of the main corrections to quenching: 
namely the stronger interquark force at small sep- 
aration (by a factor of 33/(33 — 2Nf)). How- 
ever, corrections at moderate and large separa- 
tions were treated empirically. This year, several 
groups have used dynamical quark simulations di- 
rectly. These have mainly involved 2 flavours of 
staggered quarks and the coupling has been deter- 
mined by relating ajj^ to the plaquette action by 
perturbation theory. For example [Q, the pla- 
quette action gives av(3.41/a) = 0.1785(56) at 
(3 — 5.6 for 2 flavours of staggered quarks with 
ma — 0.01. Here the error comes from assuming 
the coefficient is ±1. The scale has been set 
by either mp or the IP- IS splitting in onia. Af- 
ter processing the Nf = and Nf = 2 results to 
obtain Nf = 3, in each case the result obtained 
for 3 light fiavours at moderate q has been run 
to g = Mz in the conventional way to facilitate 
comparisons. 

Using Wilson valence quarks, the result |Q ob- 
tained is 

a^g(Mz)= 0.111(5) 
where an extrapolation of Sa to = has 
been made. The result for Nf — 2 has been cor- 
rected to Nf = 3 by using 2-loop running down 
to « 0.5 GeV with 2 fiavours and then back up 
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with three. The authors compare the quenched 
and unquenched a values and find that they over- 
lap at low q values. This is some support for their 
method of correcting to Nf = 3 (and incidentally 
for the method used to go directly from Nf — to 
Nf = 3 |29|). A similar approach using staggered 
valence quarks shows sizable flavour breaking 
effects, particularly for heavier quarks, and thus 
relatively large systematic errors. 

Another group js^ has used NRQCD to eval- 
uate the bb spectrum with dynamical quark con- 
figurations. They extrapolate in Nf from 
Nf = and 2 to get to Nf = 3. Using a fixed 
dynamical quark mass aniq — 0.01 at [3 — 5.6, 
gives 

awsi^^z) = 0.115(2) 

A preliminary result has been obtained us- 
ing the same configurations but with Wilson va- 
lence fermions. The main difference comes from 



the lattice scale as 



= 1.9 GcV rather than 2.4 



GeV as found from NRQCD. Their result is 
^Wsi^z) = 0.108(6) 

Note also that using tq to set the scale on 
the same configurations ||l^ yields ro/a = 5.2(3) 
which gives « 2.0 GeV. This agrees with 
the Wilson valence result but differs from the 
NRQCD value from the same configurations. The 
NRQCD value relies on perturbative evaluation of 
coefficients in the lagrangian, so non-perturbative 
corrections introduce an error which is hard to 
quantify. The direct full QCD determinations 
have scale errors of 0{a) which can be reduced 
by increasing f3. 

In view of the uncertainties of order 4% found 
in the quenched case with much smaller lat- 
tice spacings, it seems prudent to include errors 
at least of this size to allow for possible non- 
perturbative effects and possible effects from the 
unknown perturbative coefficients at 0{a^). Fur- 
ther, the scale error is presently some 20% which 
implies a 4% error in a{Mz)- In summary, I pro- 
pose that at least a 6% error should be ascribed 
to the lattice results. Even so, a reasonable aver- 
age of current lattice determinations would yield 
0.112(7). This is still of significance since the 
world average quoted above was 0.117(5). 

To improve on this, calculations at increased 
(3 are needed to explore non-perturbative effects 



and 0{a) effects. Further variation of ruq and 
Nf should be explored. The 0{a^) terms in the 
perturbative matching also need to be evaluated. 
Several of these steps are under way and it should 
soon be possible to evaluate at w 10 GeV 
to high accuracy from lattice studies. Assuming 
the conventional 3-loop evolution from 10 GeV, 
this will yield a[Mz) to 1%. 



4. THE INTERACTION 
COLOUR SOURCES 



BETWEEN 



As well as explicit study of the spectrum and 
coupling of hadrons, it is instructive to look at 
qualitative features of QCD. The confining mech- 
anism is the clearest example. One consequence 
of confinement is that the potential energy V{R) 
between static colour sources in the fundamental 
representation grows with increasing separation 
R. In the quenched approximation, this potential 
continues to grow approximately linearly for large 
i? - as KR where K is the string tension. The 
name string tension comes from the observation 
that a string theory can yield such a behaviour. 
To study this confining force in more detail, there 
are several topics to pursue: a detailed measure- 
ment of V{R) at large R is of interest; the spatial 
distribution of the energy in the colour flux can 
be determined; the spin dependence of the long 
range force can be explored. These topics have 
been of interest for many years. Increased com- 
putational resources have enabled major increases 
in precision recently. Because it is qualitative fea- 
tures that are to be explored, it is sufficient to 
use SU(2) colour for initial investigations. This 
allows considerable increase in precision since it 
is computationally faster. 

As well as studying the strong force between 
two heavy quarks, it is of interest to explore the 
much weaker forces between hadrons. These can 
be studied by looking at 4 quark potentials in 
SU(2). 

4.1. The self-energy of the hadronic string 

The self-energy of closed string-like flux tubes 
of length L (called torelons) has been studied in 
SU(2) gauge theory For i > 1 fm, the en- 
ergy of such torelons was found to have a compo- 
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nent j'KjiL with / = 1.00(3), exactly the coeffi- 
cient expected from string fluctuations. This is a 
test of the zero point fluctuations of a string. The 
string excitations themselves give rise to excited 
states of the potential between static sources and 
they have been found previously to agree well. 

4.2. The colour flux tube 

The Wuppertal group |Q have explored the 
spatial distribution of the colour fields between 
static quarks in SU(2) with separations out to 
i? « 2 fm. At smaller i?, their results are similar 
to features seen previously [|l|. The question of 
interest is the transverse profile of the colour fiux 
tube. They find that the width of this transverse 
flux profile increases quite fast with R until i? « 1 
fm when the increase levels off. By looking at 
the differences in the distribution of colour flux 
as i? is increased, they find consistency with the 
picture that a string-like configuration exists for 
i? > 0.75 fm. The data are shown in fig ^ and are 
consistent with the expected increase of the width 
as log i? in this region, but the evidence is not 
compelling. The breakdown into colour magnetic 
and colour electric contributions is measured and 
this is a constraint on phenomenological models 
of the confining flux tube. 

4.3. The spin dependent force 

From earlier work, we know that the force be- 
tween static quarks has a short ranged compo- 
nent roughly like one gluon exchange and a long 
range component which is dominantly in the Y{ 
rather than VJ^ spin-orbit term (this corresponds 
to Lorenz 'scalar' confinement in the usual de- 
scription). In order to get more precision, it is 
necessary to go to smaller lattice spacing to study 
the spin dependent potentials since they vary very 
rapidly at small K. 

Recent work confirms the earlier qualita- 
tive results and shows some evidence for small but 
significant departures from the simple picture of 
one gluon exchange together with a string com- 
ponent. The region studied is i? < 0.3 fm. In this 
i? region, the new results confirm that the con- 
fining component is dominantly in V[. As well as 
this string tension term K in V/(i?), there is also 
evidence for an additional attractive 1/i?^ piece. 
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Figure 9. The half width p of the energy and 
action flux tube at source separation K in units 
of the string tension K. 



Since 1 gluon exchange does not contribute to V/, 
such a component could arise from the string fluc- 
tuation term (-7r/12i? in Vb(i?)) or from multiple 
gluon exchange. 

A careful theoretical study of the nature of the 
spin-orbit force to be associated with the long 
range string fluctuation component would be very 
useful since some relevant lattice data now exists. 

4.4. The force between hadrons 

The colour force between static fundamental- 
representation quarks is far stronger than the 
residual 'strong' force between hadrons. It has al- 
ways been difficult to explore this 'strong' force in 
lattice simulation because the binding energy is so 
much smaller than that from the colour force. For 
example for SU(2) of colour and for 2 quarks and 
2 antiquarks, the difference t4 — 2V2 between the 
4-body potential and that of the two pairs taken 
separately is typically much less than 1% of 2V2. 
Nevertheless, results have been obtained for the 
binding energy for a number of different geometri- 
cal arrangements of the 4 static sources |Q . The 
value of the observed binding energy can be un- 
derstood qualitatively from mixing between dif- 
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ferent pairings of the two body potentials. One 
strong conclusion from this study is that the 
four body system is not well approximated as a 
sum of two body components. 

5. CONCLUSIONS 

Often the progress in a field is hard to see year 
by year. New ideas become accepted gradually 
and less accurate methods wither gently. The 
progress is often that problems are no longer dis- 
cussed, since new methods have circumvented the 
obstacle. 

One area of progress has been in understanding 
the best way to study the QCD coupling on a 
lattice. The agreement between different methods 
gives confidence. Present estimates correspond to 
a value ay^{Mz) = 0.112(7) and a clear route 
exists to reduce this error to about 1%. 

Quenched QCD is currently in good shape. 
The sources of error are understood and steps 
are being taken to reduce them. Accurate spec- 
tral results in the quenched approximation are of 
importance, both because they enable us to un- 
derstand QCD better and because they provide 
a calibration for calculations of matrix elements, 
etc. The full QCD program is less advanced be- 
cause of computational limitations. Studies so far 
have used rather small lattices, rather coarse lat- 
tice spacing, and rather heavy quarks. With fu- 
ture increases in computer power and future im- 
provements in algorithms, it will be possible to 
explore full QCD as thoroughly as the quenched 
approximation is currently studied. 
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